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1.  INTRODUCTION 


In  many.instances,  the  engineering  staff  is  requested  by  program  planning  personnel  to 
provide  answers  to  a  certain  design  question  regarding  a  launch  vehicle,  which  is  in  the 
development/quahfication  stage,  or  to  abnormal  behavior  of  a  vehicle  in  flight.  Often,  a 
quick  response  is  important,  in  order  that  pertinent  decisions  can  be  made  which  have  major 
impact  on  cost  and  scheduling  of  a  program.  This  kind  of  request  is  common  in  the  area 
concerning  flows  inside  solid  rocket  motors.  Existing  techniques  for  motor  internal  flow 
analysis  usually  require  time-consuming  and  extensive  labor  for  geometry  modeling  and 
flow-field  solution  and  cannot  provide  analysis  results  in  an  e^qpeditious,  sati^actory  manner. 

In  this  study,  an  efficient,  intelligent  method  for  the  solution  of  compressible  viscous  flows 
inside  solid  rocket  motors  is  developed  to  provide  timely  analysis  support  of  launch  vehicle 
design,  qualification,  and  flight  hardware  evaluation  for  the  Air  Force  Space  and  Missile 
Systems  Center.  The  method  applies  an  adaptive,  unstructured,  finite-element  mesh 
generation  technique  and  a  cell-centered,  second-order,  flux  correction  tran^ort  scheme  for 
the  solution  of  Navier-Stokes  equations,  which  include  every  term  in  the  viscous  stress  tensor. 
The  method,  originally  developed  at  the  University  of  >\Wes,  Swansea,  England  (Ref.  1),  and 
NASA  Langley  Research  Center  and  contractors  (Refs.  2  and  3)  for  the  study  of  a  plane 
two-dimensional,  supersonic  flow  over  a  circular  cylinder,  is  extended  and  streamlined  to  the 
solution  of  axisymmetric  flows  inside  solid  rocket  motors  of  arbitrary  configuration  with  a 
nonuniform  inlet  boundary  condition. 

With  known,  initial  motor  interior  configuration,  the  program  developed  in  this  study 
allows  the  computational  boundaries— including  propellant  bum-backs,  ablated  motor  case 
insulation,  and  eroded  nozzle-exit  cone  suri^ce  at  any  time  slices  after  ignition— to  be 
determined  intelligently.  The  basic  meshes  and  corresponding  compressible  viscous  flow 
fields  then  are  generated  for  the  specified  computational  boundaries.  Based  on  the  results  of 
flow  analyses  on  the  basic  meshes,  the  element  sizes  are  refined  in  high-gradient  regions  and 
unrefined  in  low-gradient  regions.  The  remeshes,  thus  obtained,  are  appropriate  for 
cost-effective  production  flow-field  runs.  The  process  of  calculation  is  fully  automated, 
resulting  in  smooth  flow  in  transmitting  computational  results  from  one  module  to  another. 

The  advantages  of  using  the  present  method  are  as  follows: 

•  Easy  implementation  of  an  exact  design  configuration 

•  Optimum  use  of  computational  resources 

•  High  resolution  analysis  results  for  minimum  cost 

•  Elimination  of  elapsed  time  through  an  automated  procedure 

•  Same  procedure  for  any  complicated  configuration 

•  Same  mesh  generation  program  for  different  applications,  such  as  thermal  and 
structural  analyses 

•  Significant  reduction  in  effort  required  for  a  complete  solid  rocket  motor  internal 
flow-field  solution 

It  can  be  said  that  because  of  its  flexibility  and  versatility  in  treating  complicated  geometry 
in  multidimensional  space,  the  adaptive,  unstructured  finite-element  method  for  the  solution 
of  Navier-Stokes  equations  will  be  the  main  computational  technique  for  flow  analysis  in  the 
coming  decade.  In  the  following  sections,  the  formulation  and  the  solution  method  will  be 
discussed,  implication  of  the  method  to  the  solution  of  axisymmetric  compressible  viscous 
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flows  inside  the  Jet  Propulsion  Laboratory  (JPL)  nozzle,  the  Titan  IV  solid  rocket  motor 
upgrade  (SRMU)  nozzle,  and  the  Star-48  motor  with  a  very  complicated  chamber  and  nozzle 
geometry  wHl  be  illustrated. 
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2.  FORMULATION 


The  Navier-Stokes  equations  in  cylindrical  coordinates  for  a  plane  two-dimensional  or  an 
axisymmetric  space  are: 


+  ip 

dx  dr  r 


utrx  +  vx„  -J-  K — 


QU 

P  +  QU^ 
guv 

u(p  +  Qe\ 


VXn  +  UXn  +  K- 


QV 

guv 

P  + 

yip  +  geX 


where  8  =  0  or  1,  depending  on  the  configuration 

f  Tn- = V-v  + 

3  dr 

2  ^  -  ,r,  du 

^xx  =  V-v  -I-  2p— 

13  dx 

1  2  (5v 

I  ^ee  =  V-v  +  2p— 


^fx  =  '^xr  =  P 


du  ^  dv  \ 
dr  dx  j 


dv  dv  du 

V-v  = - 1 - 1 - 

dr  r  dx 
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Tq  +  5i 
r  +  5i 


Sutherland’s  theory  of  viscosity  for  air  is 


where  To  =  492*R,  Si  =  199®R,  and  4o  =  2.4865  x  10"^  Ibf  -  sec/in.^. 
The  power  law  for  viscosity  of  combustion  gas  is 


i^ere  A  »  constant  for  a  combustion  gas. 

Multiplied  by  a  differential  area  (2irr)^  drdx,  and  integrated  over  an  element,  with  the 
^^lication  of  divergence  theorem  on  the  first  two  terms  of  the  right-hand  side,  Eq.  (1) 
becomes 


vdiere 

Tc 

r 

E„ 

e 


-  area  of  element 

s  radial  coordinate  of  element  centroid 
=  boundary  of  element 
=  flux  normal  to  the  boundary  of  element 
==  subscript  indicates  flow  variable  at  element  centroid 


(2) 


3.  SOLUTION  METHOD 


Notice  that  in  Eq.  (2),  8  =  0  is  for  a  plane  two-dimensional  configuration,  and  8  =  1  is  tor 
an  axisymmetric  configuration.  When  governing  equations  are  written  in  the  form  shown  in 
Eq.  (2),  the  second-order,  flux  correction  tran^ort  scheme  for  .solution  of  the  plane 
two-dimensional  Navier-Stokes  equations,  contained  in  the  computer  program  LARC/NESS 
(Langley  Adaptive  Remeshing  Code  and  Navier-Stokes  Solver)  and  developed  in  References 
1  throu^  3,  can  be  modified  easily  and  adopted  to  axisymmetric  configurations  by  including 
the  terms  required  for  axisymmetry.  The  detailed  solution  method  based  on  this  scheme  is 
discussed  in  Reference  2. 

For  rocket  nozzle  flows,  the  initial  condition  for  this  time-dependent  scheme  is  derived 
from  a  one-dimensional  isentropic  assumption.  The  exit  boundary  condition  is  based  on  a 
simple  linear  extrapolation,  since  the  flow  is  assumed  to  be  supersonic  at  the  exit  plane,  and 
the  error  generated  firom  the  extrapolation  is  not  e}q>ected  to  propagate  back  and  affect  the 
upstream  results.  There  are  several  different  schemes  (e.g..  Refs.  4  and  5)  for  implementing 
an  inlet  boundary  condition  in  a  time-dependent  nozzle  flow  calculation.  The  characteristics 
formulation  shown  in  Reference  4  and  used  in  Reference  6  provides  a  fairly  smooth  subsonic 
flow  in  the  physical  domain  and  is  applied  here  in  an  adaptive,  unstructured  mesh  form,  if  the 
inlet  boundary  is  not  on  the  propellant  burning  surface.  When  the  inlet  boundary  lies  on  the 
propellant  burning  surface,  the  flow  variables  on  the  inlet  boundary  are  determined  from 
propellant  bum-rate,  chamber  temperature,  <md  pressure. 

For  a  viscous  flow,  the  nozzle  wall  is  a  non-slip  boundary  with  a  specified  wall 
temperature.  For  an  inviscid  flow,  the  nozzle  wall  is  a  slip  boundary  with  the  velocity  vector 
tangent  to  the  wall  surface.  For  a  cell-centered  scheme,  such  as  the  one  used  in  this  study,  the 
radius  of  element  centroid  is  always  a  finite  value  and  never  becomes  zero.  Therefore,  no 
singularity  problem  is  associated  with  the  calculation  of  flow  variables  at  the  elements 
containing  the  centerline  as  one  of  the  element  boundaries  in  the  present  cell-centered 
formulation.  The  flux  across  the  centerline  boundary  is  zero.  For  other  finite-element 
methods  with  a  node-centered  scheme.  Reference  7  provides  a  technique  for  resolving  the 
numerical  difficulty  encountered  at  the  singular  centerline.  From  numerical 
e?q)erimentation,  it  is  advisable  to  use  a  smaller  integration  step  size  in  an  axisymmetric  space 
than  that  in  a  plane  two-dimensional  space. 

The  plane  two-dimensional  LARC/NESS  program  (Ref.  2)  has  been  enhanced  for  an 
axisymmetric  calculation  in  this  study  through  the  following; 

•  Developing  an  auxiliary  program  for  intelligently  determining  computational 
boundaries  that  include  propellant  bum-backs,  ablated  motor  case  insulation,  and 
eroded  nozzle-exit  cone  surface  at  any  time  slice  after  ignition 

•  Streamlining  the  program  to  automatically  blend  a  structured  mesh  in  the  boundary 
layer  with  an  unstructured  mesh  in  the  main  flow  region  for  any  arbitrary  geometry 

•  Employing  an  isentropic  initial  guess  to  start  a  nozzle  flow  calculation  on  an 
unstructured  mesh 

•  Implementing  procedures  for  evaluating  nonuniform  inlet  boundary  conditions 

•  Including  every  term  in  the  governing  Eq.  (1)  for  an  axisymmetric  configuration 

•  Incorporating  a  power  law,  viscosity-temperature  relationship  for  combustion  gas, 
in  addition  to  Sutherland’s  theory  of  viscosity  for  air 
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•  Simplifying  the  meshing  and  remeshing  procedure  through  automated  data 
transmissions  between  the  different  computational  modules 

•  Utilizing  publication-quality  plotting  routines  for  post-processing  data  displays 
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4.  THE  JPL  NOZZLE 

The  flow  inside  the  JPL  axisymmetric  convergent-divergent  nozzle  (Ref.  8)  with  a  45“ 
entrance  and  a  15  *  exit  straight  wall  tangent  to  a  circular  throat  (with  ratio  of  throat  radius  of 
curvature  to  throat  height  =  0.625)  provides  a  classic  comparison  for  transonic  nozzle  flow 
analysis.  The  cold  flow  test  was  conducted  with  air  at  a  stagnation  pressure  of  70  psia  and  a 
stagnation  temperature  of  540*R. 

Figure  1  shows  the  basic  mesh  and  remesh  for  the  JPL  nozzle  flow  calculation.  The 
remesh  was  obtained  after  500  integration  steps  on  the  basic  mesh.  There  are  20  stretched, 
structured  elements  (with  the  first  element  height  10"^  in.)  covering  the  thickness  of  the 
viscous  layer,  which  envelops  the  nozzle  wall  boundary.  The  total  number  of  elements  is  3298 
(1140  structured  quadrilateral  +  2158  unstructured  triangular  elements)  for  the  basic  mesh 
and  2328  (1280  structured  quadrilateral  -H  1048  unstructured  triangular  elements)  for  the 
remesh. 


REMESH 


Figure  1.  JPL  Nozzle  Meshes 

Sutherland’s  theory  of  viscosity  for  air  is  employed  in  the  JPL  nozzle  flow  calculation.  The 
computed  dimensionless  pressure  ratio  (P/Po)  and  Mach  number  contours  are  shown  in 
Figure  2.  The  pressure  contour  has  a  constant  increment  of  0.025,  and  the  Mach  number 
contour  has  a  constant  increment  of  0.1.  The  flow  recompression  near  the  wall  downstream 
of  the  throat,  which  was  observed  in  a  cold  flow  test  (Ref.  8),  can  be  seen  from  the  contour 
plots.  The  enlarged  view  of  the  remesh  and  Mach  number  contours  in  the  throat  region  is 
depicted  in  Figure  3,  where  the  structured  mesh  near  the  wall  blending  with  the  unstructured 
mesh  for  the  main  flow  is  distinguishable.  The  calculated  boundary  layer  thickness  shown  in 
Figure  3  is  less  than  that  observed  in  Reference  8,  mainly  because  of  the  laminar  flow 
assumption  in  the  present  study. 
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PRESSURE  CONTOUR 


MACH  NUMBER  CONTOUR 

Figure  2.  JPL  Nozzle  Pressure  and  Mach  Number 


MACH  NUMBER  CONTOUR 


Figures.  JPL  Nozzle  Throat  Region 
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The  computed  pressure  and  Mach  number  distributions  along  the  wall  and  centerline  are 
illustrated  in  Figures  4  and  5.  Shown  in  the  same  figures  for  comparison  are  nozzle  wall 
geometry  and  test  data  from  Reference  8.  The  computed  results  agree  very  well  with  test  data 
in  the  entire  subsonic-transonic-supersonic  flow  regime.  This  indicates  that  the  boundary 
conditions  are  implemented  correcdy  for  the  axisymmetric  nozzle  flow  calculation  with  an 
adaptive,  unstructured  mesh  and  assures  further  application  of  the  computer  program  to 
other  nozzle  configurations.  In  general,  the  present  viscous  analysis  results  show  similar  flow 
features  to  those  presented  in  Reference  6.  A  somewhat  better  flow  resolution  is  obtained 
from  the  present  analysis  than  that  of  the  invisdd  flow  analysis  of  Reference  6  in  the 
low-^eed  compression  comer  area,  where  the  viscous  effect  prevails. 


Figure  4.  JPL  Nozzle  Pressure  Distribution 


Figure  5.  JPL  Nozzle  Mach  Number  Distribution 
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No  test  data  are  available  for  comparison  with  the  computed  convective  heating  rate 
shown  in  Figure  6.  In  the  present  calculation,  the  nozzle  wall  temperature  is  held  constant  at 
500 *R.  The  analysis  indicates  that  gas  pressure  stays  essentially  unchanged  until  it 
approaches  the  throat.  The  highest  positive  (gas  to  wall)  heating  rate  in  Figure  6  is  indicative 
of  high  gas  energy  recovery  near  the  wall  at  a  high  contraction  ratio  location.  As  flow  passes 
through  the  throat,  the  gas  pressure  decreases  rapidly  and  the  gas  temperature  falls  below  the 
constant  value  at  the  wall,  and  a  negative  (wall  to  gas )  heating  rate  occurs.  The  absolute  value 
of  a  computed  heat  exchange  between  the  gas  and  the  nozzle  wall  for  laminar  flow  has  a  peak 
at  the  throat,  which  also  is  observed  from  the  turbulent  flow  heating  rate  measurement 
carried  out  in  Reference  9  for  a  convergent-divergent  nozzle  flow.  The  energy  equation  in 
Eq.  (2)  converges  much  slower  than  the  continuity  and  momentum  equations.  An  accurate 
heating  rate  calculation  requires  many  integration  steps  with  very  fine  structured  elements  in 
a  viscous  layer.  This  is  an  area  which  needs  to  be  improved  in  the  future. 


-3.0  -2.5  -2.0  -1.5  -1.0  -03  0.0  03  1.0  1.5  2.0 

X  (m.) 


Figure  6.  JPL  Nozzle  >^^11  Heating  Rate 
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5.  THE  TITAN  IV  SRMU 


The  Air  Force  Titan  IV  SRMU  (Ref.  10),  shown  in  Figure  7,  is  being  developed  to  launch 
large  payloads.  This  is  a  126-in.  diameter,  112-ft  long,  three-segment  motor  with  a  graphite 
epoxy  composite  case  (33,600  lb  lighter  than  a  steel  case).  The  motor  is  loaded  with 
688,850  lb  hydro)Q^lterminated  polybutadiene  (HTPB)  propellant  (69%  ammonium 
perchlorate  +  19%  aluminum  +  12%  binder)  and  weighs  about  772,750  lb.  The  nozzle 
throat  is  made  of  graphite/phenolic,  and  the  exit  cone  has  a  tape-wrapped  carbon/phenolic 
forward  insulator  and  silica/phenolic  aft  insulator.  The  nozzle  has  a  composite  elastomer 
flexseal  with  a  maximum  6”  gimbal.  The  maximum  mass  flow  rate  is  5700  Ib/sec  which 
produces  approximately  1.6  million  Ibf  thrust  for  each  SRMU  during  liftoff.  The  Titan  IV 
with  two  SRMUs  on  each  side  of  the  core  vehicle  is  designed  to  put  a  41,000  lb  nominal 
payload  into  low  earth  orbit.  There  will  be  five  full-scale  static  motor  tests  during  the 
development  and  qualification  phases. 


Figure  7.  Titan  IV  SRMU  Baseline  Design 

Ordinarily,  a  cold  flow  test  would  be  required  to  obtain  flow-field  information  for  a 
thermostructural  evaluation  of  nozzle  design  adequacy.  The  cold  flow  test  would  cost  at  least 
one-half  million  dollars  and  would  require  a  one-year  set-up  time  for  a  limited  number  of 
propellant  burn-backs.  The  efficient,  intelligent  method  developed  in  this  report  for  the 
solution  of  viscous  flows  inside  solid  rocket  motors  can  generate  flow  solutions  for  any 
number  of  propellant  bum-backs  within  a  few  days. 

The  expected  motor  head-end  chamber  pressure  history  is  given  in  Figure  8.  Figure  9 
shows  ablated  motor  case  insulation,  eroded  nozzle,  and  propellant  bum-backs  for  the  Titan 
IV  SRMU  at  0, 25, 50, 75,  and  100  sec  after  ignition.  The  motor  case  insulation  ablation  rate 
is  5  mil/sec.  The  erosion  rate  varies  along  nozzle-exit  cone  stations  and  is  projected  to  be 
10  mil/sec  at  the  throat.  The  propellant  bum-rate  is  a  function  of  pressure  (0.03454*P®  ^^). 
For  the  SRMU  flow-field  calculation,  the  ratio  of  specific  heats  is  1. 156,  the  Pfandtl  number  is 
0.57,  and  the  power  law  viscosity-tenmerature  relationship  has  an  ejq)onent  A =0.67  and 
reference  chamber  viscosity  1.382«10^1bf-sec/in.2  at  temperature  6364*R. 
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Figure  8.  SRMU  CSiamber  Pressure  Histoiy 


Figure  9.  SRMU  Nozzle,  Motor  Case,  and  Propellant 
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Figure  10  shows  the  basic  meshes  and  remeshes  at  0, 50,  and  100  sec  after  ignition.  The 
remeshesare  obtained  after  500  integration  steps  on  the  basic  meshes.  There  are  5717,6321, 
and  6733  elements  for  the  basic  meshes  at  0, 50,  and  100  sec  after  ignition.  Correspondingly, 
there  are  4442, 5838,  and  5950  elements  at  0, 50,  and  100  sec  for  the  remeshes.  These  include 
20  stretched,  structured  elements  covering  the  thickness  of  viscous  layer,  w^ich  envelops  the 
nozzle-exit  cone  wall.  Figure  11  shows  the  computed  Mach  numl^rs  and  dimensionless 
pressure  ratio  (P/Po)  contours  at  0, 50,  and  100  sec  after  ignition.  The  enlarged  view  of  the 
remeshes  and  Mach  number  contours  in  the  nozzle  region  is  depicted  in  Figure  12.  The  Mach 
number  and  pressure  distributions  on  the  wall  and  centerline  are  given  in  Figures  13  and  14, 
respectively.  Figures  15  and  16  show  the  detailed  pressure  distributions  in  the  nozzle  and  exit 
cone  regions. 

To  cover  the  entire  nozzle  and  the  exit  cone  pressure  field,  a  nonuniform  contour 
increment  is  used  in  the  Titan  IV  SRMU  pressure  contour  plots.  For  a  pressure  ratio  greater 
than  or  equal  to  0.8,  a  constant  increment  of 0.025  is  used.  For  a  pressure  ratio  less  than  0.8,  a 
decreasing  factor  of  0.9  is  applied  to  each  subsequent  contour  value.  This  covers  the  pressure 
field  better  than  that  using  a  linear  or  a  logarithmic  contour  increment.  The  Mach  number 
contour,  however,  can  be  adequately  displayed  with  a  constant  increment  of  0.1. 

The  computed  ratio  of  static-to-chamber  pressure  at  the  exit  plane  of  the  exit  cone  is  1% 
at  the  wall  and  0.6%  at  the  centerline.  For  a  viscous  analysis,  the  Mach  number  at  the  wall  is 
zero.  The  “Mw”  in  Figure  13  representing  the  gas  Mach  number  at  the  edge  of  the  boundary 
layer  is  evaluated  from  an  isentropic  relation  at  the  local  wall  pressure.  The  Mach  number  at 
the  exit  plane  of  the  exit  cone  is  3.3  at  the  wall  and  3.5  at  the  centerline.  There  are  some 
differences  in  the  flow  field  at  different  times  fi-om  ignition;  these  can  be  observed  from  the 
contour  plots  and  Mach  number  and  pressure  distributions.  It  is  interesting  to  note  that  the 
pressure  remains  nearly  constant  at  the  chamber  pressure  condition  on  the  backside  of  the 
submerged  nozzle  and  deaeases  rapidly  as  gas  flows  through  the  throat.  This  pressure 
distribution  is  important  for  detailed  nozzle  thermostructural  analysis  and  nozzle  design 
adequacy  assessment. 
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Figure  10.  SRMU  Meshes 
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Figure  11.  SRMU  Pressure  and  Mach  Number 
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Figure  12.  SRMU  Nozzle  Region 
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Figure  13.  SRMU  Mach  Number  Distribution 
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6.  THE  STAR-48  MOTOR 


The  Star-48  motor  (Ref.  1 1)  is  an  upper-stage  motor,  which  provides  the  impulse  to  propel 
a  spacecraft  (up  to  2800  lb)  from  a  low  earth  orbit  into  a  geosynchronous  transfer  orbit.  The 
motor  has  propellant  loading  up  to  4400  lb.  The  motor  case  is  made  of  titanium  and  is 
internally  insulated  with  silica-Mled  ethylene  propylene  diene  monomer  (EPDM).  The 
nozzle  is  made  of  three-directional  carbon/carbon,  and  the  exit  cone  can  have  either  a 
two-directional  carbon/carbon  or  a  carbon/phenolic  liner.  The  motor  assembly  and  the  initial 
grain  geometry  are  shown  in  Figure  17.  The  chamber  pressure  history  for  a  low-temperature 
(30 *F),  long-bum  (94  sec)  firing  is  given  in  Figure  18. 


For  this  complicated  initial  grain  and  submerged  nozzle  geometry,  the  auxiliary  program 
developed  under  this  study  intelligently  determines  the  propellant  bum-backs,  ablated  motor 
case  insulation,  and  eroded  nozzle-exit  cone  configurations  at  any  number  of  time  slices. 
These  configurations,  which  are  shown  in  Figure  19  at  0, 20, 40, 60,  and  80  sec  after  ignition, 
are  utilized  direcdy  as  the  computation^  boundaries  for  the  mesh  generation  and, 
subsequently,  for  the  flow-field  solution.  The  actual  ablation  and  erosion  data  firom  the 
post-fire  measurement  of  full-scale  static  motors  in  a  qualification  program  (Ref.  1 1)  are  used 
in  the  auxiliary  program  for  determining  the  computational  boundaries.  The  propellant 
bum-rate  is  a  fimction  of  pressure  (0.0355«P®  For  this  Star-48  motor,  the  inlet  boundaries 
lie  on  the  propellant  burning  siu^ce.s,  and  the  program  automatically  provides  the  inlet  flow 
conditions  at  any  specified  time  slice. 


Figure  20  depicts  the  basic  mesh,  remesh,  and  computed  Mach  number  contours  for  the 
initial  motor  configuration.  The  total  number  of  elements  is  6438  for  the  basic  mesh  and  3503 
for  the  remesh.  These  include  10  stretched,  structured  elements  covering  the  thickness  of  the 
viscous  layer,  wdiich  envelops  the  nozzle-exit  cone  wall.  Figure  21  gives  the  same  information 
at  80  sec  into  the  motor  firing.  The  total  number  of  elements  is  8035  for  the  basic  mesh  and 
4854  for  the  remesh.  For  the  Star-48  flow-fields  calculation,  the  ratio  of  ^ecific  heats  is 
1.150,  the  Prandtl  number  is  0.57,  and  the  power  law  viscosity-temperature  relationship  has 
an  exponent  A = 0.67  and  reference  chamber  viscosity  1.382»10^  Ibf-sec/in.^  at  temperature 
6364“R. 


The  Mach  number  and  pressure  distributions  are  given  in  Figures  22  and  23,  respectively. 
For  a  viscous  analysis,  the  Mach  number  at  the  wall  is  zero.  The  “Mw”  in  Figure  22 
representing  the  gas  Mach  number  at  the  edge  of  the  boundary  layer  is  evaluated  firom  an 
isentropic  relation  at  the  local  wall  pressure.  At  0  sec,  the  Mach  number  at  the  exit  plane  of 
the  exit  cone  is  3.91  at  the  wall  and  4.11  at  the  centerline.  The  computed  ratio  of  the 
static-to-chamber  pressure  at  the  exit  plane  of  the  exit  cone  is  0.29%  at  the  wall  and  0. 14%  at 
the  centerline.  At  80  sec,  the  Mach  number  at  the  exit  plane  of  the  exit  cone  is  3.85  at  the  wall 
and  4.01  at  the  centerline.  The  computed  ratio  of  static-to-chamber  pressure  at  the  exit  plane 
of  the  exit  cone  is  0.34%  at  the  wall  and  0.20%  at  the  centerline.  At  80  sec  after  ignition,  the 
computed  pressure  and  Mach  number  distributions  near  the  throat  are  not  smooth,  ^ich 
reflects  the  effect  of  the  unevenly  eroded  nozzle  wall  on  the  flow  behavior  in  that  region.  The 
change  in  nozzle  surface  configuration  firom  0  to  80  sec  afifects  the  Mach  number  distribution 
and  pressure  load  on  the  wall.  On  the  backside  of  the  submerged  nozzle,  the  pressure  remains 
essentially  constant  at  the  chamber  condition. 
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Figure  17.  Star-^8  Motor  Assembly 


Figure  18.  Star-48  Chamber  Pressure  History 


Figure  19.  Star-48  Motor  ConjBguration 
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Figure  20.  Star-48  Meshes  and  Mach  Number  at  0  sec 
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Figure  22.  Star-48  Mach  Number  Distribution 
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Figure  23.  Star-48  Pressure  Distribution 

The  purpose  of  the  calculation  with  the  Star-48  motor  given  here  is  to  illustrate  the 
capability  of  the  adaptive,  unstructured  finite-element  method  for  an  efficient  solution  of 
viscous  compressible  flows  inside  a  solid  rocket  motor  with  a  complicated  grain  and  nozzle 
configuration.  A  coarse  mesh  has  been  assigned  in  the  viscous  layer  to  prevent  unnecessary 
consumption  of  the  computational  resource  for  this  demonstrated  calculation.  Performing 
500  integration  steps  on  the  basic  mesh  of  6438  elements  takes  19  min  on  a  Cray  XMP/18 
supercomputer.  Performing  500  integration  steps  on  the  remesh  of  3503  elements  takes 
11  min  on  the  same  computer.  Note  that  remeshing  results  in  a  significant  reduction  in  the 
total  number  of  elements  required  to  adequately  cover  the  entire  physical  domain  and  in  the 
amount  of  computation  time.  A  meaningful  flow-field  solution  can  be  obtained  after 
approximately  2000  integration  steps  on  the  remesh.  Since  the  Star-48  is  a  spinning  rocket 
motor,  final  analysis  based  on  the  technique  developed  in  this  study  requires  a  fine  mesh  in  a 
viscous  layer  and  consideration  of  the  particle-phase  subject  to  radial  and  axial  acceleration 
in  the  flow  field. 
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7.  CONCLUSION 


An  efficient  method  has  been  developed  for  the  solution  of  axisymmetric  compressible 
viscous  flows  inside  solid  rocket  motors  of  arbitrary  geometry.  A  gas-particle,  two-phase 
viscous  flow  analysis,  similar  to  that  given  in  References  6, 12,  and  13  for  an  inviscid  flow,  will 
be  the  subject  of  study  in  the  coming  months.  Further  study  will  be  directed  toward  the 
solution  of  a  gas-partide,  two-phase  viscous  flow  with  radial  acceleration  (spinning  rocket 
motor)  and  axial  acceleration  (flight  motor)  bocty  force  terms  in  the  governing  equations. 
Adaptive,  unstructured  finite-element  methods  for  three-dimensional  flow  analyses  have 
been  reported  in  the  literature  (Re&.  14  and  15).  A  similar  technique  for  a  three-dimensional 
canted  or  gimbaled  nozzle  with  consideration  of  turbulent  flow  modeling  also  is  a  subject  of 
future  study. 
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